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Abstract. We calculate the propagator of the domain wall fermion (DWF) of 
the RBC/UKQCD collaboration with 2+1 dynamical flavors of 16^ x 32 x 16 
lattice in Coulomb gauge, by applying the conjugate gradient method. We 
find that the fluctuation of the propagator is small when the momenta are 
taken along the diagonal of the 4-dimensional lattice. Restricting momenta in 
this momentum region, which is called the cylinder cut, we compare the mass 
function and the running coupling of the quark-gluon coupling as,gi{q) with 
those of the staggerd fermion of the MILC collaboration in Landau gauge. 

In the case of DWF, the ambiguity of the phase of the wave function 
is adjusted such that the overlap of the solution of the conjugate gradient 
method and the plane wave at the source becomes real. The quark-gluon 
coupling as^g^iq) of the DWF in the region q > l.SGeV agrees with ghost- 
gluon coupling as{q) that we measured by using the configuration of the MILC 
collaboration, i.e. enhancement by a factor (1 -|- c/q"^) with c ~ 2.8GeV^ on 
the pQCD result. 

In the case of staggered fermion, in contrast to the ghost-gluon coupling 
as{q) in Landau gauge which showed infrared suppression, the quark-gluon 
coupling as^giiq) in the infrared region increases monotonically as (7 — > 0. 
Above 2GeV, the quark-gluon coupling Us^g^ {q) of staggered fermion calcu- 
lated by naive crossing becomes smaller than that of DWF, probably due to 
the complex phase of the propagator which is not connected with the low 
energy physics of the fermion taste. 



1 Introduction 

In the calculation of quark-gluon vertices in the infrared region, non-perturbative 
renormalization is possible by calculating the quark propagator in a fixed gauge. 
The calculation of the quark propagator on lattice is reviewed in [1]. In our 
previous paper [2], we studied the quark propagator of staggered fermion in 
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Landau gauge using the full QCD configurations of relatively large lattice (24^ x 
64) of MILC collaboration [4j available from the ILDG data base [3J. In the last 
year, full QCD configurations of the domain wall fermion (DWF) of medium size 
(16^ X 32 X 16) were released in the ILDG and in this year large size (24'^ x 64 x 16) 
were released [5] from the RBC/UKQCD collaboration [6]. 

In these configurations the length of the 5th dimension was fixed to be 16. 
In this paper we show the results of the medium size DWF configurations and 
compare with the results of the large size staggered fermion. 

Charcteristic features of infrared QCD are confinement and chiral symmetry 
breaking. The confinement is related to the Gribov copy i.e. non gauge unique- 
ness, which makes the sharp evaluation of physical quantities difficult, and we try 
to fix the gauge in the fundamental modular region [7]. Chiral symmetry break- 
ing is speculated to be related to instantons [8]. The Orsay group discussed that 
the infrared suppression of the triple gluon coupling is due to instantons [91, IIQJ. 
The running coupling from the quark-gluon coupling in quenched approximation 
also showed similar infrared behavior [11]. Our simulation of the ghost-gluon 
coupling in Landau gauge obtained by configurations of the MILC collaboration 
showed infrared suppression, but in Coulomb gauge the running coupling ai{q) 
of MILC and of RBC/UKQCD did not show suppression [12j. Thus, it is in- 
teresting to check the difference of the Coulomb gauge and the Landau gauge, 
staggered fermion and DWF, and the ghost-gluon coupling and the quark-gluon 
coupling. 

The domain wall fermion (DWF) was first formulated by Kaplan in 1992 
[13\ [H] by assuming that the chiral fermion couples with the gauge field in the 
fifth dimension. The model was improved by Narayanan and Neuberger [15] and 
Shamir [TBI HZ] ) such that the gauge field are strictly four dimensional and are 
copied to all slices in the fifth dimension. The model was applied in the finite 
temperature simulation of 8^ x 4 lattice with Lg from 8 to 32 lattices [18j and to 
quenched simulation of 8^ x 32, 12^ x 32, and 16^ x 32 lattices with Ls from 16 
to 64 [20j. 

The fermionic part of the Lagrangian formulated for the lattice simulation is 

[Muniisi], 

Sf{'>P,^/J,U) = - ^ '4'xADF)x,s;y,s'i'y,s', (1) 
x,s;y,s' 

where 

{DF)x,s;y,s' = <5s,s'-D|y + S^^yD^,,,. (2) 

The interaction Z)!' contains the gauge field and the interaction in the fifth di- 
mension defined by does not contain the gauge field [20j. 

The bare quark operators are defined on the wall at s = and s = — 1 as 

qx = Phi^xfl + Pr'(Px,Ls , (3) 

1 + 75 1 ~ 75 

where Pr = — - — and Pl = — - — are the projection operator. For the Dirac's 

7 matrices, we adopt the convention of ref. [18], in which 75 is diagonal. 
In the DWF theory, a Lagrangian density in the fermion sector 

= - i^)^ + ^{MPr + M^Pl)^) (4) 
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with an operator M acting on the left-handed and right-handed field was pro- 
posed [15j. In this method, the free fermion propagator becomes 

On the lattice, one introduces the bare quark mass nif that mixes the two 
chiralities, and 5 dimensional mass M5. The lattice simulation of DWF propaga- 
tor is performed by introducing a Hamiltonian, whose essential idea is given in 
Appendix. This formalism was adopted in the Schwinger model [22] and in the 
4-dimensional lattice simulation [18\ [23l [6] . 

Instead of using the transfer matrix method, we calculate the quark prop- 
agator by using the conjugate gradient method in five dimensional spaces. We 
interpret the configuration at the middle of the two domain walls in the fifth di- 
mension as the physical quark wave function. We measure the propagator of the 
domain wall fermion using the configurations of the RBC/UKQCD, and compare 
the propagator with that of the configurations of MILC [2] . We measure also the 
quark-gluon coupling from the quark propagator, by applying the Ward identity. 

The organization of this paper is as follows. In sect. 2, we present a formula- 
tion of the lattice DWF and its numerical results are shown in sect. 3. In sect. 4 
the lattice calculation of the staggered fermion propagator which we adopted in 
[2] is summarized and in sect. 5, a comparison of the DWF fermion-gluon and 
the staggered fermion-gluon is given. Conclusion and discussion are given in the 
sect. 6. Some comments on the Hamiltonian is given in the Appendix. 



2 The lattice calculation of the DWF propagator 

In this section we present the method of calculating the DWF propagator. 

Using the Dp defined in eql2]we make a hermitian operator Dh = j^R^Dp, 
where (i?5)ss' = 5s^Ls-i-s' is a reflection operator as 



Dh 



I5PR 75p" - 1) 

V 75(dII - 1) 75PL 



75PR 75iD\\ - 1) \ 
75PR 75(^11-1) I5PL 



I^Pl 



-m 



H^Pr ) 



where Pril = (1 i Ts)/^, and -1 in (Z)" — 1) originates from D^^,. 
The quark sources are sitting on the domain walls as 

q{x) = Pl^{x, 0) + Pb^{x, L, - 1). (6) 

We take Pl^{x,s) oc e~^r)^ , Pr^{x,s) oc e"(^' r ^)^, (s = 0, 1, • • • , - 1) with 
r = 0.842105, which corresponds to (r/L^)/(l + r/Ls) = 0.05 [24j- 
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In the case of quenched approximation, a condition on the Pauh-Villars reg- 
ularization mass M for producing a single fermion with the left hand chirality 
bound to s = and the right bound to s = — 1 is < M < 2. In the free the- 
ory, the condition that the transfer matrix along the 5th dimension be positive 
yields a restriction < M < 1 |15j . However, in a quenched interacting system, 
M = 1.8 was adopted [19] and since in the conjugate gradient method there is 
no M < 1 constraint, we adopt the same value. 

We define the base of the fermion as 



^{x) = \(t>L{x, 0), (I)r{x, 0), • • • , (I)l{x, Ls - 1), (I)r{x, Ls - 1)) 



(7) 



where * means the transpose of the vector, and the (pL/ni^j h) contains color 3x3 
matrix, spin 2x2 matrix and Ux >^ riy x Hz x nt site coordinates. We measure 
9 • 4 matrix elements on each site at once. 

The wave functions (j)i/ji{x, s) are solutions of the equation 
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-1 



4'l{x,s) 
-B 

-ct 



c 

-B 



(f)L{x,s) 

4'r{x,s) 



-B C 
Ct B 



(f)L{x,s) 

(I>r{x,s) 



where 



and 



Here 



I 4 



I 4 

c'= oE(^/^(^)'5-+A,y-t^V(yM- 



M5 = Me{s - Ls/2) 



-M 
M 



s < 
s > 



2 

Ls-l 



(8) 



(9) 



(10) 



:ii) 



is the mass introduced for the regularization. Using the 7 matrices as defined in 
ITSl. we obtain 



Cix,y)PR 
C\x,y)PL 



I 4 

9 '^iUf,{x)5x+f,-y - ^,{y)5x-f,,y)^PR, 
U^{y)5x-^,y)^^PL 



4 

\ ^1 



where 



2 ^ 



-i(T2 

-I J 



(12) 



and Z'^^ 



\ -I J 
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The conjugate gradient method for solving the 5 dimensional DWF propa- 
gator is a simple extension of the method we used in the staggered fermion [2j, 
since the degrees of freedom in the 5th dimension can be treated as if they are 
internal degrees of freedom on each 4 dimensional sites. 

As in the transfer matrix method, we define 

n (14) 



5-M; 



-Dueo 



and 



such that 



U=[l ^f"" ] (15) 



where even-odd decomposition is done in the 5 dimensional space. We solve the 
equation for 

Hi 

using the source 

The solution on the odd sites is calculated from that of even sites as 

In the process of conjugate gradient iteration, we search shift parameters 
for [27] for c/)i^ and aj} for and in the first 50 steps we choose ak = 
Min{a^, ) and shift 0^^^ = (t^k ~ '^k4'k <^fe+i = 4>k ~ ^k4>k ™ ^he last 
25 steps we choose ak = Max(a^,a^), so that the stable solution is selected for 
both (pL and (pR. 

The convergence condition attained in this method is about 0.5 x 10~^. One 
can improve the condition by increasing the number of iteration, but the overlap 
of the solution and the plane wave do not change significantly. 

To evaluate the propagator, we measure the trace in color and spin space of 
the inner product in the momentum space between the plane waves 

x{p) = \xl{p,o),xr{p,o), ■ ■ ■ ,xl{p,Ls - 1),xr{p,Ls - 1)) 
and the solution of the conjugate gradient method 

!f = \<pLip, 0), Mp, 0), • • • , - 1), Mp, Ls - 1)) 
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as 

Tr{xip,s)PL^ip,s)) = Zb{p){2N,)Bl{p,s), 

Tr(x(p, s)Pr^{p, s)) = Zb{p){2N,)Br{p, s) (19) 

and 

Ti{x{p,s)ii>PL^{p,s)) = ZA{p)/{2Nc)ivAL{p,s), 

T¥(x(p, s)i^PR^ip, s)) = ZAip)/{2N,)ipABip, s) (20) 

1 2TTpi 

where Pi = - sm — — [pi = 0, 1, 2, • • • , Ni/2). 
a Ni 

On the lattice at each s the 4-dimensional torus is residing. We perform the 
Fourier transform in the 4-dimensional space, but take the momentum in the 5th 
direction to be zero since it corresponds to the lowest energy state. Z^ip) and 
Zb{p) are the wave function renormalization factor. 

When p4 = 0, the term B{p, s) are given by the matrix elements of {xR-, ^l) 
and {xl,^r)- The operator ^ yields matrix elements of {x,^^l) and {Xi^'^r)- 
The propagator is parametrized as 

0(^^-1 -i^ + M\p) -i^ + Mjp) 

^P^ + M{p)MHp) ^^ ^p^ + M^p)M{p) ^ ' 



where 



and 



^^^^_Re[Bn{p,Lj2 



M\p) 



Re[AR{p, Ls/2)] 

Re[BL{p,Ls/2)] 
Re[AL{p,L,/2)y 



2 TTpi 

The momentum assignment Pi = — sin is introduced for removing doublers 

a Ni 

using the Wilson prescription. 

The Ai{p) has zero eigenfunction and dim{KerM) = n/j = 1 and the M\p) 
does not have zero eigenfunction and dim{KerM'^) = = 0. 



3 Numerical results of the DWF propagator 

The configurations of the RBC/UKQCD collaboration are first Landau gauge 
fixed and then Coulomb gauge fixed {diAi = 0) as follows [12], We adopt 

the minimizing function Fu[g] = HA^lp = Ex,i (^3^,*^^^,^) , and solve 

d^Ai{x, t) = using the Newton method. We obtain e = — qjj'^^^^ from the eq. 

diAi + diDi{A)e = 0. Putting ^(a;, t) = in U9i{x,t) = g{x,t)Ui{x,t)g^ {x + i,t) 
we set the ending condition of the gauge fixing as the maximum of the diver- 
gence of the gauge field over A'^^ — 1 color and the volume is less than 10~^, 
MaXx^a{diAx^i)"' < 10~^. This condition yields in most samples 
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We leave the remnant gauge on Aq{x) unfixed, but since the Landau gauge 
preconditioning is done, it is not completely random. We leave the problem of 
whether a random gauge transformation, or the remnant gauge fixing on Ao{x) 
modify the propagators, but we do not expect drastic corrections will happen. 

Using the gauge configurations of RBC/UKQCD collaboration after Coulomb 
gauge fixing, we calculate Tr(x(p, s)0l(p, s)) and Ti{xip, s)i;^(pL{p, s)) and 
Tr(x(P! s)4>r{p, s)) and Ti{x{p, s)i;^4'R{p, s)) at each 5-dimensional slice s. Num- 
ber of samples is 49 for each mass ruf = 0.01, 0.02 and 0.03. We measured in 
certain momentum directions of nif = 0.01, 149 samples. 

In our Lagrangian there is a freedom of choosing global chiral angle in the 
5th direction, 

^ ^ e'^T^^tP, ij; 'ilje-''^^tp. (22) 

We adjust this phase of the matrix element such that both Tr{x{p, 0)(j)L{p, 0)) 
and Tr{x{p, Lg — 1)^r{p, Lg — 1)) are close to a real number. Namely, we define 



e 



i0L 



T>(x(p,O)0l(p,O)) 
\Tr{x{p,0)MP,m' 

-iOn ^ T^{x{p,Lg-l)MP,Ls-l)) 
\Tr{x{p,Ls-l)MP,Ls-m 

and sample-wise calculate e*^ such that 

is minimum. When p is even and the momentum is not along the diagonal of the 
four dimensional system, we also calculate e^^ such that 



9l^W _ i|2 + |e^^«e-'^' - Ip (24) 



is minimum, but the final results by multiplying e*'' and e*'' are similar. 

In the calculation oi B^/n, we define matrix elements multiplied by the phase 

as ^.—^ 

{x{p,s)<Pl{p,s)) = {x{P,s)(l)L{p,s))e~''^, 

{x{P, s)(I)r{p, s)) = {x{p,s)4>R{p,s))e'-''. 

and correspondingly denote Bl/ji{p,s) multiplied by the phase e~*^ and e*^ as 
Bl/r{p,s), respectively. 

In the calculation of Ai/Rip, s), we diagonalize 

[{xiPx, s)4'LiPx, s))ai + (xiPy, s)(pL{Py, s))a2 

+{x{Pz, s)(Pl(.Pz, s))<73 + {x{Pz, s)(I)l{Pz, s))il] (25) 

and 

[{X{PX, s)(t)R{px, S))ai + {x{Py, s)(t)R{py, S))(72 

+{x(.Pz, s)4>r{Pz, s))a3 + {xiPz, s)(t>RiPz, s))il], (26) 
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where I is the 2x2 diagonal matrix, and define the A^iJiip, s) multiphed by the 
phase as s). 

The term B^jfi is a sum of color-spin diagonal scalar, while the term A^jr 
is a color-diagonal but momentum dependent spinor and we take the positive 
eigenvalue. 

In order to minimize the artefact due to violation of rotational symmetry 
of the lattice we restrict the momentum configuration to be diagonal in the 4-d 
lattice. This prescription which is called cylinder cut [28] is already adopted in 
ghost propagator [3Dj and in quark propagator [21 [T] calculations. 

In general, there is a mixing between (j)L and (/>/j and there is a sign problem 
i.e. the sign of Re[BiiJi{p,Ls/2)\ and i?e[^x//j(p, Ls/2)] becomes random. The 
sign is related to the sign of the source at s = and s = — 1. But, when the 
cylinder cut is chosen, the sign problem does not seem to occur. 

In the calculation of the propagator of DWF, the mass originates not only 
from the mid-point matrix Q^"^P^ defined as 

Q^'^^K,s' = Pl5s,Ls/2^s',Ls/2 + Pr^s,Ls/2-i5s',Ls/2-1 (27) 

but also from Q^^'^ defined as 

Q^^\^,, = PiSsfiSs'fi + PrSs,Ls-iSs'l.-i (28) 

At zero momentum the numerator Bl{p = 0,s = 0) becomes 1 and it gives 
a contribution of mjQ^^^ = ruf. Since there is no pole mass in (/)/j(s,Zs), the 
value of Br{p = 0, s = — 1) is not physical. In the midpoint contribution 

Re[BL/R{p, Ls/2)] 

= , we take into account that the numerator of the mass function 

Re[AL,a{p,Ls/2)] 

contains {2Nc) x {2Nc) coherent contributions and divide by the multiplicity. 

The FiglUis the mass function of nif = 0.01/a = O.OlTGeV. The momenta 
correspond to p = (0, 0, 0, 0), (1, 1, 1, 2), (2, 2, 2, 4), (3, 3, 3, 6) and (4, 4, 4, 8). The 
dotted lines are the phenomenological fit 

. , cyl2"+^ ruf , , 

Since the pole mass Q^'*"^ is not included in the plots, ruf is set to be here. The 
corresponding values of 0.02/a = 0.034GeV and 0.03/a = 0.050GeV are similar. 

In the fit) choose a equals 1,1.25 and 1.5 and searched best values for 
c and A. We found the global fit is best for a = 1.25. The fitted parameters are 
given in Tabled! 

In the case of staggered fermion in Landau gauge, we adopt in this work 
the staple plus Naik action on mj = 0.0136GeV and 0.027GeV configurations [21 
[29] . We fixed the parameter a = 1 and obtained A = 0.82GeV and 0.89GeV, 
respectively. In general A becomes larger for larger a, but A of RBC/UKQCD 
seems larger than that of MILC, which is also observed in the quenched overlap 
fermion propagator [3T]. In the case of MILC, A becomes smaller for smaller mass 
ruf, but in the case of RBC/UKQCD, it is opposite. Analytical expression of the 
quark propagator in Dyson-Schwinger equation is formulated in [32] [33] and a 
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•mud/ a 


rus/a 


c 


yl(GeV) 


a 


DWFoi 


0.01 


0.04 


0.24 


1.53(3) 


1.25 


DWF02 


0.02 


0.04 


0.24 


1.61(5) 


1.25 


DWF03 


0.03 


0.04 


0.30 


1.32(4) 


1.25 


MILC/i 
MILC/2 


0.006 
0.012 


0.031 
0.031 


0.45 
0.43 


0.82(2) 
0.89(2) 


1.00 
1.00 



Table 1. The fitted parameters of mass function of DWF(RBC/UKQCD) and staggered 
fermion (MILC) witli tiie staple plus Naik action. 

comparison with these lattice data are given in |341 [35] . The mass function of the 
staggered fermion is close to that of the Dyson-Schwinger equation of Nf = 3, 
but larger than that of the A'^^ = 0. 

0.6 



0.5 - 




p(GeV) 

Figure 1. The mass function in GeV of the domain wall fermion as a function of the modulus of 
Euclidean four momentum p(GeV). nif = 0.01. (149 samples). Blue disks are m_t (left handed 
quark) and red boxes are mn (right handed quark). 

The error bars are taken from the Bootstrap method after 5000 re-samplings 
|36l I37j . The re-sampling method reduces the error bar by about a factor of 10 
as compared to the standard deviation of the bare samples. 

We measured also momentum points (p, 0,0, 0), (0,p, 0,0), (0, 0,p, 0) and 
{p,p,p,0) with p = 1,2,3 and 4, but the error bars of the mass function are 
found to be large especially at (2,2,2,0). We observed systematic difference of the 
magnitude of mass functions of {p, 0, 0, 0) with even p and odd p. Such problems 
would be solved by improving statistics, and systematically correcting deviation 
from the spherical symmetry, but at the moment these problems can be evaded 
by adopting the cylinder cut. 

4 The lattice calculation of the staggered fermion propagator 



The lattice results of fermion propagator using staggered (Asqtad) action and 
overlap action are reviewed in [Ij. We calculated the propagator of staggered 
fermion of the MILC collaboration using the conjugate gradient method [2]. The 
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0.6r 
0.5 
0.4 



0.3 
0.2 
0.1 





V 



2 3 

p (GeV) 



Figure 2. The mass function in GeV of the staggered fermion MILC/2 as a function of the 
modulus of Euclidean four momentum p(GeV). The staple plus Naik action is adopted 



inverse quark propagator is expressed as 



'ti)al3 



■ sin(pp 



24 



sin(3pp 



+ md, 



(30) 



where a„ = 0, 1, 13 u 



staggered fermion takes values 

27rm,, 



0,1 and dap = D/, '^ap/3^|mo<i2- The momentum of the 
+ vra^ where 

L, 



m,, 



0, 



The 7 matrix of staggered fermions is 

(7/.)a/3 = (-1)' 

where 



(m) 



if u < fi 
otherwise 



(31) 



(32) 



(33) 



The A{p) is defined as 



i^^i-ir-p^Tr[J2Sa(3{p)] = l6N,p^A{p) (34) 

where 16 is the number of taste. 

The staggered fermion incorporating the lattice symmetries including parity 
and charge conjugation by introducing a general mass matrix is formulated in 
|42j . In our model, we do not incorporate the general mass matrix, but take the 
same mass as MILC collaboration. 

The chiral symmetry of the staggered fermion of the MILC collaboration is 
currently under discussion [441 [l5l I46j . In our simple model, the charge conju- 
gation operator can be taken as C = —^74 and 6*7^(7^ = —j/i- We interpret p 
in the original as q in crossed channel and since staggered actions are invariant 
under translation of 2a, we modify the scale by a factor of 1/2: 



1 



9 1 

-sin(p^) - — sin(3p^) 



(35) 
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5 A comparison of the DWF-gluon and the staggered fermion-gluon 
coupHng 

In this section we calculate the running coupling in the crossed channel qq — > 
gluon of the DWF and compare with that of the staggered fermion. In Coulomb 
gauge, the quark gluon vertex from the three point Green function becomes 
iQlil] 

G^.{p,q) = ! <fx ! d'ye'Py-'^P+'i^^myMOhMOM^))- (36) 



When the momentum transfer q is small, the vertex function satisfying the Ward 

d 

identity Zyr^{p) = —i— — S^^{p) becomes 

OPfj, 

r^{p,q) = S-\p)G^{p,q)S-\p) 

= ^"''[giip'^H. + i92{p^)q^, + 9^{p^)p^,i] 

(37) 

where 5°'^ is the delta function in color space and Zy is the vertex renormalization 
factor. 

When the contribution of the ghost-quark coupling [25] is ignored, the vector 
current Ward identity allows us to extract the running coupling as,g^{q) from 
the difference of 5"i(p + f ) and S^Hp - f ) l26j 

- i[S-H{p + f ),|0) - S-'iip - |),|0)] = Z^'Aoip)^. (38) 

When the crossing is performed, the momentum transfer becomes p + ^ — 

{—p + ^) = 2p. For massless fermion with cylinder cut, Ap^ = p^ and thus p can 
be interpreted as q. 

In the case of DWF, we diagonalize 

Y^iiAL'^^ip + |), - ^L°^(P - |),>,] (39) 

and ^ 

Y,[{Ar^^{p + |), - Ar^P{p - |),>,], (40) 
i=i 

and to get the running coupling, we evaluate the average and multiply the nor- 

malization Z^ (p) oc -z— — = — „f , x — where „, , is from normalization of 
E{p) 2E{p) 2 2E{p) 

the B = A4A in the original direct channel, which should be proportional to 

-— and — comes from fixing the incoming wave as o or o and Z = (2N..)'^ 

2E{p) 2 ^ B y y \ c) 

comes from the relative normalization of A and B. 
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In Figini we show the running couphng of DWFqi and MILC/i. An enhance- 
ment of Us^g^iq) of RBC/UKQCD above 2GeV region could be the effect of the 

condensate due to instantons [TO]. Using operator product expansion, the 
Orsay group fitted the lattice data above 2.6Gev as 



a 



Latt 



(41) 



where the parameter c is proportional to the condensate. They obtained 



c = 2.7(1.2) ^^J^^lMg^pii^GeV 

We show pQCD result with N j = 3 without (dot-dashed line) and with the 
A? condensate effect (dashed line) where c = 2.8GeV^ is used |17J as in the 
analysis of the ghost-gluon coupling. Consequence of the condensates on the 
mass gap and quark confinement is discussed in [48]. 

In g > 2GeV, the running coupling as,g^{(l) of MILCji is smaller than that 
of RBC/UKQCD. We think it is due to the complex phase of the staggered 
fermion in the ultraviolet region which is not related to the low energy physics 
of the fermion taste [l6]. In the previous analysis of propagator of the MILC/2 
in which the bare s-quark mass is close to the bare ujd quark mass, we observed 
an anomalous behavior when the Asqtad action is adopted [29]. Since the sample 
size was not large, we cannot exclude the possibility that the anomalous behavior 
of the staggered fermion disappear in the simulation of larger number of samples. 
We leave these problems in the future. 

The running coupling in the infrared region is consistent with the experimen- 
tal data extracted by the JLab group [l3] . They compared the proton form factor 
and the neutron form factor and by adopting the Drell-Hearn-Gerasimov sum 
rule in infrared and the Bjorken sum rule in ultraviolet, extracted the running 
coupling. 



0.1560) 
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Figure 3. The running coupling as,gi(g) of MILC/i (blue disks) and DWFoi (green points). 
The dash-dotted line is the pQCD result and the dashed line is the pQCD with the A? conden- 
sate contribution. The red points are data extracted from the experiment by the JLab group. 
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6 Conclusion and discussion 

The mass function and the running couphng in Coulomb gauge of the gauge 
configuration of RBC/UKQCD collaboration were calculated and compared with 
those of the staggered fermion. We adopted the conjugate gradient method and 
imposed a reality condition on the overlap of the distorted wave and the plane 
wave at the position of the fermion sources. 

We observed that the quark-gluon coupling a^.gi {q) of the DWF is consis- 
tent with the ghost-gluon coupling as{q) of MILCji in q > l.SGeV region. The 
running coupling as,gi{q) of the staggered fermion in Landau gauge (MILC/i) 
does not show infrared suppression, in contrast to the ghost-gluon coupling [37]. 
The discrepancy of the ghost-gluon coupling as{q) and the quark-gluon coupling 
'^s,gi{q) in Landau gauge suggests that there is a problem in the ghost-gluon 
coupling, and/or the color structure of the loop given by a product of ghosts 
[381 [39] . The difference of Landau gauge quark-gluon vertex in quenched config- 
uration [11] and in unquenched configuration that we measured, suggests impor- 
tant contribution of fermions in the dynamics of gluons. Orsay group interpreted 
the infrared suppression of triple gluon vertex of quenched configuration as the 
instanton effect [41] . In the infrared, however, the vacuum amplitude with pres- 
ence of instantons contain zero-mode divergence [8J. In the expression of the 
triple gluon vertex 



it was assumed that G^^\q^ , q^ , q^) oc -^(P^^ilPf) ^^'^ G'^^^(g^) = 

— {p^I{qp)^), n being the instanton density and p being the instanton radius. 
8 

When p is large, the zero-mode divergence could overwhelm q'^ dependence and 
yields a constant as{q). 

In a supersymmetric theory, the zero-mode divergence from fermion and from 
boson are shown to cancel out [l9]. In quenched simulation of gluonic systems, 
the fermionic zero mode divergence is absent, and consequently incorrect large p 
dependence of instantons could have introduced the infrared suppression of the 
triple gluon running coupling. 

Although we do not consider the supersymmetric Yang-Mills theory, we show 
in Appendix that M and could be regarded as supersymmetric interactions. 
Phenomenologically, the cancellation of zero-mode divergence from quark field 
and that from gluon field in the conjugate gradient calculation of the quark 
propagator seem to have introduced the correct infrared behavior. 

The new method of deriving the quark propagator is encouraging, but it is 
necessary to extend the calculation to larger lattice for getting the continuum 
limit, and to extend the simulation for other momenta that are far from the 
4-dimentional diagonal axis. The origin of the fluctuation of the data outside the 
cylinder cut region is under investigation. 

Acknowledgement. The author thanks Reinhard Alkofer for a discussion on the quark propa- 
gator in Coulomb gauge and the support of author's stay in Graz in March 2008, and Hideo 
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Appendix A: The Hamiltonian of the Domain Wall Fermion 

The Hamiltonian of the free DWF can be expressed as 



and its square becomes 







^ 




(S 1) 






I) 



^^^^-^ -{i> + ^f + MHi )■ (^-2) 

Taking the eigenstates of hamiltonian including the gauge potential ^ as the expansion bases 
and identifying 

/n ^\ /nn\ 

as the supersymmetry operators that satisfy 

Q' = Qt' = 0, {Q,Q''} = H 

and [H, Q] = 0, we regard M and are a pair of supersymmetric interactions [151 150) . 

In the free fermionic theory, the number of massless right-handed particle is 
dim{Ker{M)) = nn and the massless left-handed particle is dim{Ker{M^)) = n_L. @It was 
shown that by choosing a proper operator M, one can define U = M\MM'')-^^^, such that 
U^fU ^1 and f/f/1" = M''{MM'')-^M = 1 - Q where Q is the projector on the zero eigenspace 
of M. 
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